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Abstract 

The problem of explaining the acceleration of the expansion of the universe 
and the observational and theoretical difficulties associated with dark matter 
and dark energy are discussed. The possibility that GR does not correctly de- 
scribe the large-scale structure of the universe is considered and an alternative 
gravity theory is proposed as a possible resolution to the problems. 
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1 Introduction 

The recent surprising observational discovery that the expansion of the universe is 
accelerating has led to an increasing theoretical effort to understand this phe- 
nomenon. The attempt to interpret the data by postulating a non-zero positive 
cosmological constant is not satisfactory, because it is confronted by the two serious 
issues of why the estimates from the standard model and quantum field theory lead 
to preposterously large values of the cosmological constant, and the coincidence of 
matter and dark energy dominance today 

If we simply postulate a repulsive force in the universe associated with a charge 
density, then we might expect that this force could be responsible for generating 
the acceleration of the universe. However, for a homogeneous and isotropic universe 
the net charge density would be zero, although for a finite range force with a small 
mass there will exist a non-zero charge density The effect of a Maxwell-type 
force would be to lower or raise the total energy, leaving the form of the Friedmann 
equation unchanged. Thus, we would still have to invoke exotic forms of energy with 
an equation of state, p = wp, where w is negative and violates the positive energy 
theorems. For a non-zero cosmological constant w = —1. 

In addition to the dark energy problem, we are still confronted with the puzzle of 
dark matter. Any observational detection of a dark matter candidate has eluded us 
and the fits to galaxy halos using dark matter models are based on several parameters 
depending on the size of the galaxy being fitted. The dark matter model predictions 
disagree with observable properties of galaxies 0] . 

Another problem with the dark energy hypothesis is the serious challenge to 
present day particle physics and string theory from the existence of cosmological 
horizons, which arise in an eternally accelerating universe Q. Several resolutions 



of this problem have been proposed but a cosmological horizon does produce a 
potentially serious crisis for modern particle physics and string theory. 

Challenging experimental results are often the precursors of a shifting of scientific 
paradigms. We must now entertain the prospect that the discovery of the mechanism 
driving the acceleration of the universe can profoundly change our description of the 
universe. 

Given the uneasy tension existing between observational evidence for the accel- 
eration of the universe and the mystery of what constitutes dark matter and dark 
energy, we are tempted to reconsider the question of whether Einstein's gravity the- 
ory (GR) is correct for the large scale structure of the universe. It agrees well with 
local solar system experimental tests and for the data obtained for observations of 
the binary pulsar PSR 1913+16. However, this does not preclude the possibility 
of a breakdown of the conventional Einstein equations for the large-scale structure 
of the universe. The standard GR cosmological model agrees well with the abun- 
dances of light elements from big bang nucleosynthesis (BBN), and the evolution 
of the spectrum of primordial density fluctuations, yielding the observed spectrum 
of temperature anisotropics in the cosmic microwave background (CMB). Also, the 
age of the universe and the power spectrum of large-scale structure agree reasonably 
well with the standard cosmological model. However, it could be that additional 
repulsive gravitational effects from an alternative gravity theory could agree with all 
of the results in the early universe and yet lead to significant effects in the present 
universe accounting for its acceleration 0. 

When contemplating alternative gravity theories, one is impressed with the math- 
ematical and physical robustness of GR. It is not easy to change the structure of GR 
without running into consistency problems. A fundamental change in the predictions 
of the observational data will presumably only come about from a non-trivial alter- 
ation of the mathematical and geometrical formalism that constitutes GR. From the 
cosmological standpoint, such theories as Jordan-Brans-Dicke theories of gravity 
will not radically change the Friedmann equation in the present universe. Recent 
developments in brane-bulk cosmological models have led to alterations of the 
Friedmann equation but only for the very early universe corresponding to high en- 
ergies. 

This tempts us to return to a physically non-trivial extension of GR called the 
nonsymmetric gravity theory (NGT). This theory was extensively studied over a 
period of years, and a version of the theory was discovered that was free of several 
possible inconsistencies such as ghost poles, tachyons, exotic asymptotic behaviour 
and other instabilities [|1^, |ll|, |1^, |1^ . Further research is needed to fully understand 
such problems as its Cauchy development and the deeper meaning of the basic gauge 
symmetries underlying the theory. 

As we shall see in the following, NGT can describe the current data on the 
accelerating universe and the dark matter halos of galaxies, gravitational lensing 
and cluster behaviour, as well as the standard results such as BBN, the solar system 
tests and the binary pulsar PSR 1913+16, without invoking the need for dominant. 
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exotic dark matter and dark energy. 

2 NGT Action and Field Equations 

The nonsymmetric g^u and F^^, are defined by [l^, |Tl|, [12|, |13|, |14 



and 

T _ 

The contravariant tensor g'^'^ is defined in terms of the equation 



— rf ^ + Tr^j^i. (2) 



g'^'g.. = g^'gua = 5^. (3) 

The Lagrangian density is given by 

C-NGT = C + Cm, (4) 

where 

C = g^^R.AW) - 2A^ - ^fih^'giu,] 

-igA''^iy^l^^ + gMV[^0,], (5) 

and Cm is the matter Lagrangian density (G = c = 1): 

Cm = -Sirg^-'T^,. (6) 

Here, g'^'^ = ^/—gg^^, g = Det{gfj_^), A is the cosmological constant and Rfj,u{W) is 
the NGT contracted curvature tensor: 

R,AW) = <,,^ - + Wj^p^^) - WiW^p + Wj^pW^,, (7) 

defined in terms of the unconstrained nonsymmetric connection 

w^A = rj. - ls\w,, (8) 

where 

Eq.(|D leads to the result 

= ^[f_t\] = 0. 



The contracted tensor Rij,i,{W) can be written as 

where 

The term in Eq.(|^): 

S^'Jl.U (9) 

contains the Lagrange multipher fields 0^ and a source vector J^. 
A variation of the action 

S = d'^xCT 



-NGT 

yields the field equations in the presence of matter sources: 

G^^iW) + Ag^, + S^, = 87r(r^, + K^,), (10) 

s^'''\u = -Is^'-'-^w., (11) 

Here, we have Gf^u = R^iu — ^g^iuR, and 

S,u = \f^'{9[H + Ig.J'^'^gipa] + g^^'^g^.gp.) - ^iw^w, - \g,ug''^wj¥p). (13) 

Moreover, the contribution from the variation of @) with respect to g^'^ and \/—g 
is given by 

K,. = -^M^] - Ig.Ag^'^^^JM]. (14) 

The variation of 0^ yields the constraint equations 

g''^"lj. = 0. (15) 

We have not varied the source vector J^. 
If we use ([T5|), then (|1^ becomes 

= —^Jl^^<P'^]■ (16) 

If we specify to be = (0,0,0, Jq); then (|T3|) corresponds to the three 
constraint equations 

g['°l = 0. (17) 



After eliminating the Lagrange multiplier field 0^ from the field equations (|T0D, 



we get 



where g/^*^"^^ is the Levi-Civita symbol. 
The generalized Bianchi identities 

[g"^G,.(r) + g'^"G.p(r)],, + g^^\^G,. = 0, 

give rise to the matter response equations 



(18) 
(19) 



(20) 



(21) 



A study of the linear approximation has proved that the present version of NGT 
described above does not possess any ghost poles or tachyons in the linear ap- 
proximation [O. This cures the inconsistencies discovered by Damour, Deser and 



McCarthy in an earlier version of NGT |jT5[. Moreover, the instability discovered 



by Clayton [|I^ for both massless and massive NGT in a Hamiltonian formalism, 
associated with three of the six possible propagating degrees of freedom in the skew 
symmetric sector, is eliminated from the theory. This is implemented in the NGT 
action by the covariant constraint equations ([T5|) . 



3 Cosmological Solutions 



For the case of a spherically symmetric field, the canonical form of g^i, in NGT is 
given by 
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We have 





f22) 



2^^1/2 



For a comoving coordinate system, we obtain for the velocity vector m'^ which satisfies 
the normalization condition 9{^v)U^u'^ = 1: 



u 



1 

7^' 



u = u 



0. 



(23) 



From Eq.(|T^, we get = and only the g[22,] component of g[^yT^ is different from 

= --^s^pgl'-^l (24) 



zero. The vector can be determined from 



where Sfj_ag^°"^'' = 5^- For the spherically symmetric field with = 0, it follows from 
(0) and dH) that PF^ = 0. 

The energy-momentum tensor for a fluid is 

T'"' = {p + p)u^u'' - pg^" + ^['^"1 , (25) 

where B^^'^^ is a skew symmetric source tensor. 

We can prove from a Killing vector analysis that for a homogeneous and isotropic 
universe massless NGT requires that f{r,t) = |]T^. It follows that for massless 
NGT all strictly homogeneous and isotropic solutions of NGT cosmology reduce to 
the Friedmann-Robertson- Walker (FRW) solutions of GR. For the case of massive 
NGT, it will be possible to obtain strictly homogeneous and isotropic solutions. In 
the following, we shall solve the fleld equations for a spherically symmetric inho- 
mogenous universe and then approximate the solution by assuming that the inho- 
mogeneities are small. We expand the metric g^^u) as 

g{^.u) = gfju) + ^9{,.u), (26) 

where g^J^^ denotes the homogeneous and isotropic solution of g^^u) and Sg(^f^i^) are 
small quantities which break the maximally symmetric solution with constant Rie- 
mannian curvature. We shall simplify our calculations by assuming that the density 
p and the pressure p only depend on the time t. Moreover, we assume that the mass 
parameter fi ^ and we neglect any effects due to the antisymmetric source tensor 

It is assumed that a solution can be found by a separation of variables 

a(r, t) = h{r)R\t), (3{r, t) = r^S\t). (27) 
From the fleld equations, we get 

R S \ „ , , 

where R = dR/dt and Z is given by 

_ $'p 5$(3'p a[3'p 2$ff ff 3/7 

off 2lfi 

Let us assume that Z ^ 0, then from (pHj) we flnd that R ^ S and the metric 
line-element takes the form 



ds'' = de - R\t) 



h{r)dP + P{d9'^ + sin^ 



(30) 
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From the conservation laws, we get 

If we assume that P ^ f, then the equations of motion become [jl3|, |E 

2b{r) + R{t)R{t) + 2R^{t) - R\t)W{r, t) = 47rR\t)[p{t) - p{t)], 

1 Att 
- Rmit) + -R\t)Y{t) = -R'itMt) + Mt)], 

where 

2Hr) - 



The functions W and Y are given by 



2a2/?3 2ap3 2a(3^ 2a/32 

«7/' //" 4//'/3' ^ 3/'2 



2a2^2 ct/52 a/?3 2a/32' 
/3/2 5/j2/2 3/2 4/j// // 

/33 2/54 2/32 ^ /53 /52 ■ 

Within our approximation scheme, W and Y can be expressed in the form 
_ h'f 2/2 ^ 2i?2/2 ^ 10/2 i?// 



//" 8//' ^ 3/'2 



^ 2{R' + RR)P _ lOR'p _ _3p_ 8Rff ff 



Ehminating by adding and (^Bf), we get 



R^ + h=^-^pR' + QR\ 



Ji = -Yfl(p + 3p) + iflJ-'. 
We can write Eq. (|38D as 



where 



From (|33|) we obtain 



where H = R/R, 

n = nM + nQ, (42) 

and 

^M = ^, nQ = §,. (43) 



If 6 = 0, then we get Q = 1 and 



H' = ^p + Q. (44) 



The hne element now takes the approximate form of a flat, homogeneous and 
isotropic FRW universe 

ds^ = de - R\t)[dr^ + r\de^ + sin^ ed(t)% (45) 



4 Accelerating Expansion of the Universe 

It follows from (^) that i? > when Y > Att{p + 3p). If we assume that there is a 
solution for Q and Y, such that they are small and constant in the early universe, 
then we will retain the good agreement of GR with the BBN era with prad oc 
and pm oc 1/R^. As the universe expands beyond the BBN era at the temperatures, 
T ~ 1 MeV-60 kev, then Q begins to increase and reaches a slowly varying value 
with Qq ~ 0.7 and Q%j ~ 0.3, where Q%j and Qq denote the present values of Qm 
and Qq, respectively. These values can fit the combined supernovae, cluster and 
CMB data 0. 

We observe from (^) and (^7^ that the dependence of Q and Y as the universe 
expands is a function of the behaviour of R and / and their derivatives. If / grows 
sufficiently with Rast increases, then Q and Y can dominate the matter contribution 
Pm as the universe evolves towards the current epoch. A detailed solution of the 
field equations is required to determine the dynamical behaviour of R, f, Q and Y. 
However, we can obtain some knowledge of the qualitative behaviour of Y and Q 
by assuming that / ~ /' ~ and h = 1. Then, from ( p6|) and (|37|) we obtain 

(46) 

and 

/12 • -A 

<3«S%(7. + ^«'^-^«). (47) 

If, as the universe expands, the behaviour of / is / ~ i?" with a > 2, then we can 
satisfy Y>0,Q>0, R>0 and R will increase as we approach the present epoch 
corresponding to an accelerating expansion. 

We can explain the evolution of Hubble expansion acceleration within NGT, 
without violating the positive energy conditions. Both p and p remain positive 



8 



throughout the evolution of the universe. There is no need for a dark energy and a 
cosmo logical constant. Thus, we avoid having to explain the unnatural and mysteri- 
ous "coincidence" of matter and dark energy domination. The Q contribution to the 
expansion of the universe increases at a slow rate up to a constant value today with 
Y > Anp (p ~ 0), and can then decrease to zero as the universe continues to expand, 
avoiding an eternally accelerating universe. During this evolution, the cosmological 
constant A = 0. It is then possible to avoid the existence of a cosmological horizon 
and the problems it produces for quantum field theory and string theory 



5 NGT and Dark Matter 

Galaxy dynamics observations continue to pose a problem for gravitational theories 
and cosmology. The data for spiral galaxies are in sharp contradiction with Newto- 
nian dynamics, for virtually all spiral galaxies have rotational velocity curves which 
tend towards a constant value. The standard assumption is that dark matter exists 
in massive, almost spherical halos surrounding galaxies. The standard hypothesis 
is that about 90% of the mass is in the form of dark matter and dark energy and 
this explains the flat rotational velocity curves of galaxies. This explanation is not 
economical, for it requires three or more parameters to describe different kinds of 
galactic systems and no satisfactory model of galactic halos exists [Q. 

A possible explanation of the galactic rotational velocity curves problem has been 
obtained in NGT [0. A derivation of the motion of test particles yields the total 



radial acceleration experienced by a test particle in a static spherically symmetric 



gravitational field for r > 0.2 kpc, due to a point source (we reinsert G and c) |T9 



,, GM ACc^ exp(-r/ro) / r, ,^ 



where A, C and tq = l/fi are constants which remain to be fixed. 
We choose C oc \/M and fix A to give 



where Goo is defined to be the gravitational constant at infinity 



Goo = Go(l + ^^) (50) 

and Go is Newton's gravitational constant. 

These formulas were applied to explain the flatness of rotation curves of galaxies, 
as well as the TuUy-Fisher law ||20|, L ~ v^, where v is the rotational velocity of a 
galaxy and L is the luminosity. A derivation of v gives 
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where Lq = 3Mq. For distances less than 0.5 — 4 kpc, the standard Newtonian law 
of gravity will apply. For ro = 25 kpc and Lq = 250 x 10^^ Lq, an excellent fit to 
spiral galaxies was found [|1^, 0. Moreover, a good fit to the TuUy-Fisher law was 
also obtained. 

Consider now the giant spiral galaxy M31 and our Galaxy in the local group. The 
center of M31 is approaching the center of the Galaxy at a velocity ~ 119 km/sec. 
The total mass of the local group should be very large, namely, the mass-to-light 

ratio should be ~ 100 ^M/L©^ . This big ratio is normally explained using the dark 

matter hypothesis. The distance between M31 and the Galaxy is ~ 700 kpc, so the 
additional exponential force is vanishingly small, but what is left is the renormalized 
gravitational constant. Thus, the gravitational acceleration becomes 

«(^) = — (52) 

where M* ~ 17M and from the observed mass-to-light ratio: 

V J ' ^^^^ 

we predict 







M fM\ 
— 



This agrees with the estimated ratio for luminous matter without using the dark 
matter assumption. 

Gravitational lensing effects can also be accounted for in NGT. We find for the 
angle of defiection A0, obtained in the post-Newtonian approximation 



4.G„(i + ^JMoIm\ 



A* = -^^a (55) 

where M is the mass of the galaxy and R is the distance between the galaxy center 
and the defiected ray. This prediction is close to the one obtained from the dark 
matter hypothesis. 

If we calculate the acceleration expected in our solar system, we obtain 

^^^a-aNGT^lv^/^y_ 

flNewton 2 M \rQ J 

For solar and terrestrial experiments, we find for Tq ~ 25 kpc, 5a < 10^^'^, which is 
too small a deviation from Newton's law to be detected with current experiments. 
We must still acount for the estimated value of VLm ~ 0.33 ±0.053 [^. Measure- 



ments of total baryon density give Vtsh? = 0.020 ± 0.001 [Q. For h ~ 0.7 ,we get 
~ 0.04. The usual hypothesis states that cold dark matter particles contribute 
^CDM ~ 0.3. We do expect that there is some dark matter in the universe in the 
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form of dark baryons and neutrinos with non- vanishing mass ( < 10% ). It remains 
to be seen whether an alternative gravity theory such as NGT can provide an ex- 
planation for the discrepancy between visible baryon matter and dark matter. At 
the era of structure formation, the contribution of Qq could be of order Qq ~ 0.3, 
growing to its present day value of Qq ~ 0.95 so that ~ 1. NGT would then 
be required to explain the formation of galaxy structure without CDM. These are 
issues that require further investigation. 
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